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ETUDE ECOLOGIQUE DES COLLEMBOLES MUSCICOLES 
DU SIDOBRE (TARN). 


П, Modèle mathématique de la distribution 
des espéces sur un rocher 


раг L. Bonnet, P. CassaGnau et D.C. de IZARRA. 


L'étude des Arthropodes édaphiques n’a cessé de se développer 
depuis les vingt dernières années. L’impulsion donnée à la Biocé- 
notique des Collemboles par le regretté D” GisiN dès 1943 a joué 
à cet égard un grand rôle dans l'orientation écologique des recher- 
ches en Europe. 

L'analyse des facteurs de répartition des espèces ou des biocé- 
noses s’est dégagée de façon progressive de l’aspect essentiellement 
intuitif des premiers travaux, pour gagner en précision grâce à 
l’utilisation de méthodes quantitatives permettant de donner à 
l’écologie du sol des bases plus solides. L'introduction par l’un de 
nous dès 1958 de méthodes stochastiques dans l’analyse du peuple- 
ment des Thécamoebiens du sol, et les résultats obtenus [ BONNET, 
1964] permettaient d’espérer que la pédozoologie pourrait un jour 
tirer grand profit de la généralisation de telles techniques, ceci 
dans le cadre plus vaste d’une biologie quantitative désormais en 
plein essor. 

Nos prévisions semblent pleinement se confirmer aujourd’hui 
où l’outil statistique figure dans l’arsenal méthodologique de maint 
laboratoire d'écologie édaphique. 

L'étude stochastique globale de la répartition des espèces dans 
l’espace en fonction des facteurs de l’environnement ne semble 
cependant pas avoir été entreprise jusqu’à présent en France chez 
les Aptérygotes du sol, les méthodes statistiques étant jusqu'ici 
réservées à l’analyse de problèmes restreints, ou simplement testées 
dans des mises au point de méthodologie pure. C’est la raison 
pour laquelle nous avons entrepris depuis quelques années l’étude 
des peuplements d’Aptérygotes en divers points du Sud-Ouest de 
la France et dans les pays méditerranéens. Les données ainsi 
accumulées nous permettent aujourd’hui d'envisager la répartition 
écologique de ces formes sur des bases statistiques qui viennent 
confirmer certains points encore imprécis d’études géographique- 
ment limitées réalisées précédemment par l’un de nous [Cassa- 
GNAU, 1961, 1965]. 


128 COLLEMBOLES DU SIDOBRE (2) 


Le travail ci-dessous est un premier exemple de cette tentative 
de quantification et de réification des facteurs de peuplement à 
partir des données brutes. 

Si l’engouement observé de nos jours pour la mathématisation 
de disciplines primitivement intuitives révéle un désir évident de 
se protéger contre certains pièges tendus par l'imagination, il faut 
se garder cependant de laisser prendre à l'outil un développement 
tel qu’il en arrive à masquer les problèmes qu'il prétendait primi- 
tivement résoudre; d’autre part, pour aussi élaborées que soient 
les méthodes statistiques, encore faut-il qu’elles soient utilisées 
sur des données écologiquement valables, si le but poursuivi, bien 
entendu, reste l’approfondissement des problèmes écologiques, et 
non plus l’utilisation à tout prix de techniques « dans le vent ». 

C’est la raison pour laquelle nous pensons que l’écologiste doit 
en définitive rester seul juge dans le choix des problèmes et l’ex- 
ploitation des hypothèses. 


INTRODUCTION. 


Le travail qui va suivre a pour but de montrer le rôle que peu- 
vent jouer des micro-conditions locales sur la distribution des 
espèces de Collemboles dans un espace très restreint, en l’occur- 
rence la surface d’un rocher en sous-bois: en effet la localisation 
spatiale des mousses saxicoles reflète bien souvent des conditions 
très variées de drainage et d’évaporation des eaux de pluie, ainsi 
que l’indépendance relative des peuplements vis à vis des strates 
édaphiques proprement dites (litières, humus...). 

L'intérêt d’un tel biotope réside essentiellement dans le rôle que 
va jouer l’humidité comme facteur limitant. Il est évident que ce 
rôle sera d’autant plus facile à mettre en évidence que l’on étudiera 
ce facteur dans des conditions où ses valeurs successives risquent 
de constituer des barrières écologiques graduelles pour les orga- 
nismes, et non plus dans des milieux où ces variations ne recoupent 
jamais les limites vitales des espèces considérées (horizons à satu- 
ration constante, milieu cavernicole par exemple pour la plupart 
des micro-arthropodes). 

A cet égard, les peuplements muscicoles sur rochers, dans le 
Sidobre, offrent une grande variété de niches écologiques différentes 
[IZARRA, 1969]. 

Résultant de la décomposition en boule du granite à la surface 
du plateau, le rocher étudié ci-dessous est situé en sous-bois clair 
de chënes, hétres et pins au niveau du Lac du Merle (Tarn) à 
600 m d'altitude environ, à une centaine de mètres à l'Est de la 
route départementale n° 30. 
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Fic. 1. — Schématisation du rocher et localisation des relevés. 


Profondément ancré dans le sol, il a la forme d’un demi-ellip- 
soide deux fois plus long que large, dont le grand axe serait orienté 
Sud-Nord; la zone de contact avec le sol forme une demi-ellipse 
de 3 m de petit axe sur 3 m pour le demi-grand axe (fig. 1). 

La hauteur du rocher est de 2 m environ au-dessus du sol en 
son point le plus élevé. 

Nous y distinguerons une face Nord formée par une ligne de 
faite arrondie subhorizontale au sommet et descendant en pente 
douce vers le Nord (30° environ). De part et d’autre de cette 
ligne (a), les flancs Est et Ouest ont une pente de 45° chacun (c 
et d). La face Sud est pratiquement verticale (b). 

La quasi totalité du rocher est recouverte par un tapis de mousses 
de 2 4 5 cm d’épaisseur reposant sur une trés faible couche de 
matiére organique. C’est cette mousse qui a été prélevée le 11 no- 
vembre 1967 en divers points. Les prélèvements (de 150 à 200 стз) 
sont effectués A divers niveaux sur les pentes du rocher. (1 a 7, 
pente a; 10 à 13, pente b; 15 à 18, pente c; 20, pente d). Les 
prélèvements 7, 13, 18 sont situés dans une zone déjà en contact 
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avec les strates édaphiques superficielles et peuvent donc d'une 
part bénéficier d’une remontée capillaire d’eau en période de sèche- 
resse, d’autre part être colonisés à partir des horizons du sol. 

En outre 4 prélèvements ont été faits dans les couches édaphiques 
superficielles (8, 14, 19) ou profondes (9). 

Les températures notées au moment du prélèvement étaient les 
suivantes : 


— air ambiant du sous-bois : 7,5°, 
—— mousses sur rocher : 8°, 

— humus au pied du rocher : 9,5, 
— sol profond : 10°. 
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Fic. 2. — Teneur en eau des 20 relevés au 11 novembre 1967. Blane : mousses 


sur rocher: hachures obliques : mousses en contact avec le sol; hachures 
verticales : litière - humus; pointillé : sol profond. 
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La teneur en eau des 20 prélèvements a été mesurée en labo- 
ratoire à l’aide d’une balance d’hygrométrie à rayons infrarouges. 
La figure 2 résume les valeurs trouvées pour chaque prélèvement. 
On notera en particulier la diminution progressive de la quantité 
d’eau retenue au fur et à mesure de l’augmentation de la pente du 
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Fic. 3. — Détail de la composition faunistique des différents relevés. 
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26 espèces de Collemboles ont été récoltées dans ces biotopes. 


Ce sont : 


Hypogastruridae : 
Ceratophysella armata, 


Xenylla tullbergi, 
Xenylla xaviert. 


Neanuridae : 


Neanura monticola, 
Neanura muscorum. 


Isotomidae : 


Anurophorus laricis, 
Folsomia quadrioculata, 
Folsomia spinosa, 
Isotomiella minor, 
Pseudisotoma sensibilis 
Uzelia setifera, 
Vertagopus arborea. 


Onychiuridae : 


Onychiurus armatus, 
Tullbergia callipygos, 
Tullbergia krausbaueri. 


Entomobryidae : 


Entomobrya albocincta, 
Entomobrya nivalis, 
Heteromurus major, 
Lepidocyrtus curvicollis, 
Lepidocyrtus lanuginosus, 
Lepidocyrtus violaceus, 
Pseudosinella alba. 


Symphypleones : 
Dicyrtomina ornata, 
Megalothorax minimus, 
Sminthurinus aureus, 


Sphaeridia pumilis. 


La figure 3 donne le détail de leur abondance dans chaque pré- 
lévement. Les figurés sont les mémes que ceux utilisés précédem- 
ment par l’un de nous [CASSAGNAU, 1961]. Les 20 espèces retenues 
pour l'analyse statistique ont été numérotées de 1 à 20 dans l’ordre 
alphabétique global (tableau IV). 


I. ETUDES DES RELEVES FAUNISTIQUES. 


1° Mesure de l’affinité faunistique. 


Parmi les nombreux indices qui permettent de mesurer l’affinité 
faunistique des relevés, c’est-à-dire leur tendance à comprendre 
les mêmes espèces, nous avons retenu la racine carrée ® du carré 
moyen de contingence. 


Ce choix se trouve motivé par les raisons suivantes : 

— le coefficient Ф présente l’avantage de tenir compte, en sus 
des espéces présentes, de celles qui font défaut parmi les espéces 
dont la présence est possible. 

— Il a la valeur d’un coefficient de corrélation linéaire a varia- 
bles dichotomisées, c’est-a-dire ne pouvant prendre que deux valeurs, 
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correspondant à la présence ou a l’absence des espèces. La perte 
d’information qui résulte de cette dernière condition se trouve 
largement compensée par l’utilisation possible de Ф dans des ana- 
lyses ultérieures, ce que ne permettent pas les autres indices, a 
exception du coefficient de corrélation de BRAVAIS-PEARSON, d’ail- 
leurs inutilisable ici, car son emploi nécessite la connaissance des 
fréquences absolues des individus. 


Le principe du calcul de Ф, est le suivant [ Воммет, 1964 et 1969]. 


Considérons deux relevés X et Y, comportant respectivement A 
et B espèces. On a observé C espèces communes aux deux relevés. 
L'ensemble des relevés comprend N espèces au total. 

Si le nombre de relevés est suffisamment important et si les 
prélèvements ont été pratiqués de façon judicieuse, N peut-être 
considéré comme une bonne estimation du nombre total d'espèces 
possibles dans les limites de l’ensemble écologique étudié. 


Faisons l’hypothèse d’une distribution des espèces au hasard. 
Dans ce cas, le nombre théorique d’espéces communes aux deux 
relevés X et Y est égal à : 

AB 


Nous pouvons comparer à ce nombre théorique le nombre ob- 
servé C, en formant l’écart-réduit : 


C-P 
Е RSS IE ESE OT qi 
А (NTA) В (N—B) 
N3 
puis le coefficient Ф : 
a 
Se iy — 
N 


Lorsque le nombre observé C d’espéces communes est égal au 
nombre théorique, le coefficient Ф est nul. Les deux relevés peu- 
vent alors étre qualifiés de faunistiquement indifférents. Dans les 
autres cas, ® varie entre —1 et +1, et exprime alors soit une 
incompatibilité faunistique entre X et Y, soit au contraire une 
affinité entre les deux relevés. 


2. Principes de l'analyse de l’affinité faunistique. 


Nous avons fait remarquer, dans les lignes précédentes, que 
Ф correspond à un coefficient de corrélation. 

Une affinité faunistique élevée, donc une corrélation forte entre 
deux relevés X et Y indique simplement, si nous nous placons sur 
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le plan descriptif, l’existence d’une liaison entre la présence et 
l’absence de chaque espèce dans chacun des relevés respectivement. 

Or, il est logique d’imputer cette liaison à l’action de certaines 
influences s’exerçant au sein des relevés, influences que nous qua- 
lifierons désormais de Facteurs. 

Désignons par Z,, les différents Facteurs susceptibles d’affecter 
un relevé X. Ce dernier pourra être représenté par le modèle 
mathématique suivant : 

Re Be + Es + .... + Z, 

Mais ce modèle est imparfait car il ne rend pas compte de lim- 
portance avec laquelle chaque Facteur Z, intervient dans le relevé 
X, sature le relevé X. Il nécessaire d’affecter chaque Z, d'un 
coefficient mesurant la corrélation qui existe entre la variable- 
relevé X et le Facteur (coefficient de saturation). 

Le modèle prend alors la forme de l’équation factorielle suivante : 


X = aZ, + bZ, + WZ, 
avec, pour une variable-relevé donnée : 
CFB =... Ир ЕЧ 
et, pour l’ensemble des variables-relevés : 
X а?, = A, valeur propre de la matrice d’affinité. 


A partir des valeurs propres, on calcule les pourcentages de 
variance totale respectivement extraits par chaque Facteur, ce 
qui permet de hiérarchiser ces derniers. 


Les n Facteurs peuvent étre figurés par n axes orthogonaux 
gradués en saturations. Dans l’espace ainsi défini, à chaque relevé 
correspond un vecteur dont les coordonnées sont ses saturations 
respectives par les différents Facteurs. Les extrémités de ces vec- 
teurs-relevés se situent dans un hyperellipsoide, ou nuage a n 
dimensions, à l’intérieur duquel les points correspondants sont 
d’autant plus proches que les relevés qu’ils représentent sont plus 
affines. Ainsi peuvent apparaitre des « noyaux d’affinité », composés 
des relevés les plus affines faunistiquement. 

L’examen attentif de la configuration du nuage montre que les 
points-relevés se disposent suivant un certain nombre de directions 
privilégiées qui correspondent soit aux axes des Facteurs soit a de 
nouveaux axes obtenus par rotation des axes des Facteurs autour 
de leur origine. 

Il est du plus haut intérét, pour la connaissance du déterminisme 
des affinités, de rechercher la signification de ces directions ou 
« gradients » et des Facteurs qui leur correspondent. 

Cette réification peut s’effectuer par comparaison de l’ordre dans 
lequel se disposent les relevés, d’une part dans l’échelle des satu- 
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rations (pour un Facteur donné) d’autre part dans celle des compo- 
santes écologiques ou biologiques en notre possession. 

Il arrive souvent que l’exploration du milieu soit trop incomplète 
pour que l’on puisse assimiler un Facteur du modèle à une compo- 
sante réelle du milieu. Le Facteur non réifié n’en constitue pas 
moins une hypothése de travail, pouvant guider le chercheur dans 
des investigations ultérieures. 

L’analyse factorielle de l’affinité faunistique permet donc, à par- 
tir de la comparaison de simples listes d’espéces, de former un 
modèle des influences diverses qui s’exercent dans le milieu et qui 
rendent compte de l’affinité ou de l’incompatibilité spécifiques des 
différents relevés. 


3. Résultats. 


On trouvera dans le tableau I la matrice d’affinité faunistique 
des 20 relevés comparés deux a deux (données de la figure 3), et 
dans le tableau II, l’analyse factorielle de ces affinités. 

Le calcul des valeurs propres et des vecteurs propres de la 
matrice a été effectué selon la méthode de Hotelling. Nous ne don- 
nons que les 4 premiers Facteurs, qui extraient 80,2 % de la 
variance totale. 


A. — Facteur Z.. 


Les valeurs positives fortes de ce premier Facteur caractérisent 
les relevés effectués en milieu muscinal (1, 2, 3, 4, 5, 6, 10, 11, 12, 
15, 17, 18, 20). Toujours positivement saturés en Z,, mais nette- 
ment séparés du groupe précédent par une saturation plus faible, 
les relevés 7 et 13 correspondent à des Mousses en contact avec le 
sol. 

A l'opposé, dans des valeurs très faiblement positives ou dans des 
valeurs négatives de Z, se groupent les relevés édaphiques stricts 
(8, 9, 14, 19). 

Nous assimilerons donc le Facteur Z, à un gradient rendant glo- 
balement compte du caractère édaphique ou muscinal du milieu. 


B. — Facteur Z.. 


Par suite de notre connaissance insuffisamment complète des 
conditions écologiques. nous n’avons pu assigner à ce Facteur 
aucune signification cohérente. 


C. — Facteur Z.. 


On remarquera que les Mousses du flanc Nord du rocher sont 
affectées de saturations négatives en Z, alors que celles du flanc 
Sud se caractérisent par des saturations positives fortes. Cepen- 
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( : COMPOSANTE ) 
RELEVE i = s ` ) 
( : Z, : 2, : 2, : 2, ) 
( : : ; : : ) 
( : : : : ) 
( 1 : 0,785 : = 0,202 + =- 0,403 : 0 179 ) 
( 2 : 0,785 : = 0,202 : = 0,403 : 0,179 ) 
( 3 : 0,834 : 0,026 : - 0,341 : 0,042 ) 
( 4 : 0,948 : 0,038 : - 0,103 : ~ 0,084 ) 
( 5 : 0,934 2 = 0,109 : =- 0,245 : 0112) 
( 6 : 0,948 : 0,038 : = 0,103 : - 0,084 ) 
( Fi : O,310 : 0,834 : 0,018 : 0,044 ) 
( 8 £ = 0,079 3 0,760 : - 0,312 : = 0,228 ) 
( 9 : = 0,340 : 0,426 : =- 0,470 : - 0,485 ) 
( 10 : 0,771 т 0,080 : 0,311 : - 0,414 ) 
( 11 : 0,610 : = 0,063 : 0,577 : =< 0,292 ) 
( 12 : 0,489 : 0,084 + 0,786 + = 0,043 ) 
( 13 : 0,261 : 0,315 : 0,486 : 0,573 ) 
' ( 14 с = 0,409 z 0,499 : 0,082 : 0,417 ) 
( 15 £ 0,934 t = 0,109 + = 0,295 p 0,314 9 
( 16 : 0,822 : = 0,130 : 0,205 : - 0,189 ) 
( ‘7 : 0,664 : 0,211 z: — 0,079 : 0,281 ) 
( 18 : 0,864 : - 0,025 : 0,042 : 0,013 ) 
( 19 : 0,155 3 0,869 : - 0,082 : 0,034 ) 
( 20 : 0,737 š 0,384 : 0,233 : - 0,096 ) 
( : : : : ) 
io, : : : : ) 
( Valeur propre $ 9,559 3 2,893 : 2,303 : 1,279 ) 
( : : : : 
% de variance : 47,796 : 14,466 : 11,517 : 6,395 ) 
: t : : ) 
TABLEAU П. — Analyse factorielle de l'affinité faunistique des relevés. 


dant, les biolopes correspondant aux relevés 18, 16, et 20, bien 
que situés sur le flane Nord, présentent des saturations positives. 

Cette anomalie s’explique aisément si l’on considère la teneur en 
eau de ces Mousses, inférieure a 50 %, indice de conditions hydri- 
ques proches de celles que nous rencontrons sur le flanc Sud. 

La même gradation se retrouve dans les Mousses en contact avec 
le sol : le relevé 7, effectué sur le flanc Nord, n’est que faiblement 
saturé en Жу, contrairement au relevé 13 provenant du flanc Sud. 

Enfin, un phénomène du même ordre apparaît en milieu édaphi- 
que. Bien que négativement saturés, les relevés correspondants se 
différencient nettement suivant leur situation par rapport au 
rocher : les prélèvements du côté Nord (9, 8, 19) sont affectés de 
saturations faibles, alors que le prélèvement 14 effectué au bas du 
flanc Sud, se montre plus fortement saturé en Z.. 

En résumé, Z, nous apparaît comme lié à la teneur en eau et à 
l'exposition; il rend compte de l’humidité résultante du biotope. 
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Dans ce systéme, trois grands ensembles apparaissent nettement : 
les Mousses du rocher, les Mousses en contact avec le sol et le sol. 
Chacun de ces groupes se subdivise à son tour de façon nette sui- 
vant son orientation. 


Cette distribution devient encore plus claire lorsque Гоп fait 
effectuer au système OZ,, OZ, une rotation de l’ordre de 30°. L’axe 
OZ’, ainsi obtenu (fig. 4) se réifie beaucoup plus facilement que 
OZ, en gradient édaphique-muscinal. D’autre part, le coefficient 
de corrélation de rang entre les saturations en Z’, et la teneur en 
eau des biotopes (fig. 2) présente Ja valeur trés hautement signifi- 
cative de p = 0,92. П nous est done permis d’assimiler cet axe Z’, 
aun gradient d’humidité relative résultante. 


II. DISTRIBUTION DES ESPÈCES. 


1. Mesure de l'affinité cénotique. 


Pour les raisons exposées à propos des relevés, nous avons retenu 
ici encore la racine carrée du carré moyen de contingence. Le cal- 
cul de ® s’effectue de la méme facon que pour les relevés; il suffit 
d’intervertir, dans les formules que nous avons données, les relevés 
et les espèces. 


La matrice du tableau III groupe les différents coefficients d’affi- 
nité, calculés en comparant deux à deux les 20 espèces les plus 
fréquemment représentées dans les relevés. 


2. Analyse de l’affinité cénotique. 


La méthode suivie est identique a celle que nous avons utilisée 
pour les relevés. A chaque espéce correspond une équation facto- 
rielle dont on pourra obtenir les quatre premiers termes a partir 
du tableau IV. 


3. Résultats. 
A. — Facteur Z.. 


Il sépare très nettement les espèces en deux ensembles corres- 
pondant au peuplement édaphique (saturations positives) et au 
peuplement muscicole (saturations négatives), les positions extré- 
mes étant occupées par : 


Xenylla tullbergi (19), caractéristique des Mousses aérées saxi- 
coles et corticoles sur le versant atlantique de la Montagne Noire, 
d’une part; 
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( : COMPOSANTE ) 
( ESPECE e ) 
: : : : i 
: Z š 7, $ z g Z 
[ : 1 ЖЕ. “ t 3 : k ; 
: : : : 
( 1 A. laricis t = 0,545: - 0,066: 0,095: 0,739) 
( 2 D. ornata : 0,761: 0,122: 0,406: 0,401) 
( 3 E. albocincta : — 0,042; 0,755: 0,491: = 0,113) 
( 4 Е. nivalis t = 0,259: 0,452: 0,569: 0,411) 
( 5 Е. quadrioculata : 0,929: ~ 0,172: 0,072: - 0,016) 
( 6 F. spinosa : 0,872: -= 0,247: 0,105: -= 0,072) 
C 7 H. major : 0,714: 0,285: 0,419: - 0,275) 
( 8 I. minor : 0,929: -= 0,172: 0,072: - 0,016) 
( 9 L. eurvicollis : 0,656: 0,591: 0,131: 0,167) 
( 10 L. lanuginosus : 0,668: 0,534: 0,191: - 0,015) 
( 11 L. violaceus : 0,218: 0,813: 0,305: 0,083) 
( 12 M. minimus : 0,727: 0,063: 0,240: 0,405) 
( 13 N. monticola : 0,102: 0,638: 0,127: - 0,020) 
(14 P. alba : 0,923: 0,063: 0,044: = 0,259) 
( 15 P. sensibilis t = 0,360: 0,523: 0,648: 0,267) 
( 16 T. callipygas : 0,648: - 0,502: 0,408: 0,187) 
( 17 T. krausbaueri : 0,648: - 0,502: 0,408: 0,187) 
{ 18 V. arborea £ = 0,3517 = 0,106: 0,272: - 0,303) 
( 19 X. tullbergi : = 0,678: - 0,402: 0,348: - 0,027) 
( 20 X. xavieri : = 0,538: 0,393: 0,041: - 0,154) 
( : : : : S 
( t t : : 
i Valeur propre : 8,098: 3,810: 2,107: age 
: : : : 
[ % de variance t 40,489: 19,052: 10,532: 0 
: : : : 


TaBLEAU IV. — Analyse factorielle de l'affinité cénotique. 


Pseudosinella alba (14), Folsomia quadrioculata (5) et Isoto- 
miella minor (8), des strates édaphiques sur les deux versants de 
la Montagne Noire d’autre part. 


B. — Facteur Z.. 


Au sein de chacun des deux ensembles précédents, le Facteur Z, 
opère une gradation dans l’aptitude des espèces à résister a la des- 
sication du milieu, aptitude plus faible chez les édaphiques que 
chez les muscicoles. 

L’ensemble des formes édaphiques se subdivise ainsi en un grou- 
pement a Poduromorphes et Isotomidae euédaphiques (Z, négatif), 
en un groupement à Symphypléones, enfin, en un troisième grou- 
pement constitué par des Entomobryidae humicoles à écailles (Z, 
positif). 

Les formes muscicoles se répartissent de la même façon en 
Poduromorphes et Isotomidae (Z, négatif ou faiblement positif) 
et en Entomobryidae (Z, fortement positif). 

Il est donc permis d’assimiler le Facteur Z, à un gradient tra- 
duisant l'aptitude des espèces à se libérer de leur milieu pour 
gagner des biotopes plus sujets à dessication. 
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Ce Facteur est a rapprocher du Facteur Z, mis en évidence lors 
de l’analyse de l’affinité entre les relevés. Tous deux sont liés aux 
aspects hydriques du milieu. 
Une rotation du système de référence autour de son origine 
permet d'obtenir une figuration meilleure, qui соггороге pleine- 


ment l'interprétation que nous venons de donner des Facteurs Z 
et Z, (fig. 5). 
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> eee | Fa Entomobr. 
t ‘es sr muscicoles 


Pod. + Isot.' \ 
muscicoles 5 
19 i 
\ 
\ 
П \ 
0 05 Z, 
Fic. 5. — Les vecteurs-espëces dans le système des 4 Facteurs Z,, Z,, Z,, Z,. 


A chaque espèce correspond, dans le système Z,, Z, un point a otis icone 
on mène les deux axes rectangulaires O'Z, et O'Z.. “Dans ce nouveau systéme 
de référence, Геѕрёсе est figurée par un vecteur d’origine О” et dont les 
coordonnées sont les saturations еп Z, et en Z, L’échelle des saturations en 
G et Z, est donnée dans le quadrant inférieur gauche. 
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CONCLUSIONS 


Un premier point se dégage de cette analyse: la netteté de la 
réponse et l'aspect schématique des résultats, ceci malgré un échan- 
tillonnage réduit a 20 prélévements et la perte d’information décou- 
lant de la méthode utilisée (on se rappellera que les coefficients 
d’affinité ont été calculés sur le critère des présences et des absen- 
ces et non a partir des fréquences). 

On peut en conclure que les Collemboles sont d’excellents indi- 
cateurs écologiques, particulièrement sensibles aux fluctuations du 
milieu physique dans lequel ils vivent. L’alignement des relevés 
sur la figure 4 est à cet égard très démonstratif : 18-19, 7-8-9, 12- 
13-14 par exemple suivant Гахе OZ’, ou 10-11-12 suivant Гахе 
OZ’,. Notons aussi que cette schématisation est d'autant plus nette 
que les différences écologiques entre relevés sont plus tranchées : 
la face Sud du rocher, subverticale, montre des alignements plus 
rigoureux que la face Nord, au niveau des mousses en particulier. 

La séparation entre les deux peuplements, édaphique strict et 
muscicole, est bien visible sur la figure 5 et correspond à la cou- 
pure observée dans la Montagne Noire entre les synusies 1, 2, 4, 5 
et les synusies 7, 8, 9, la mousse en contact avec le sol n’ayant 
pas ici un peuplement original mais représentant plutôt une zone 
d'interférence entre les deux peuplements précédents. 

Au niveau auto-écologique, le parallélisme qui peut s’établir entre 
la dualité « versant atlantique - versant méditerranéen » dans la 
Montagne Noire et « microversant Nord - microversant Sud » sur ce 
rocher est assez frappant. Nous retrouvons un méme déterminisme 
de peuplement pour les espèces de Xenylla : Véloignement de X. 
tullbergi (19) par rapport à X. харіегі (20) suivant Гахе OZ’, 
(fig. 5) souligne l’impossibilité pour la première espèce de coloniser 
le flanc Sud, contrairement à la seconde. Les coefficients de fré- 
quences de ces deux espèces pour la Montagne Noire (CASSAGNAU, 
1965) étaient : 

X. tullbergi : versant atlantique = 20; versant méditerranéen = 0. 

X. харіегі : versant atlantique = 5; versant méditerranéen = 15. 

La localisation d’Entomobrya albocincta (3) et @E. nivalis (4) 
dans les mousses bien drainées du flanc Sud, vertical, du rocher 
s'explique par la similitude qui existe entre ce biotope et le biotope 
typique de ces espèces dans la Montagne Noire, les mousses sur 
tronc d’arbre, limitées au versant atlantique. 

Quant à la position intermédiaire d’Anurophorus laricis (1) elle 
illustre bien la mobilité de cette espèce, capable de migrations, 
localisée dans les mousses sur arbres et sous les écorces en sous- 
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bois frais, mais devenant saxicole en paysage découvert, et méme 
humicole dans les régions septentrionales. 


La figure 5 montre de part et d’autre des axes OZ, et OZ, (ou 
mieux OZ’, et OZ’,) le groupement des espéces suivant des noyaux 
correspondant a la résultante des facteurs tels qu’ils ont été définis 
plus haut : gradient muscicole — édaphique (Z’,), possibilité de se 
libérer de ces deux types de biotopes (Z’,). 

Ce n’est pas un hasard si, en fait, les affinités écologiques obser- 
vées correspondent la plupart du temps à des affinités systémati- 
ques. Ceci rejoint l'interprétation écologique des grandes tendan- 
ces phylétiques des Collemboles développée récemment par l’un de 
nous |CASSAGNAU, 1970] : nous observons d’une part une tendance 
a un mode de vie interstitiel dans les familles ayant gardé une 
métamérisation homonome, une furca parfois réduite, voire absente 
(Poduromorphes et Isotomidae), d'autre part un mode de vie libéré 
du milieu interstitiel (vie épigée ou cavernicole) dans les lignées 
à furca et appendices antérieurs bien développés, à métamérisation 
hétéronome, à recouvrement tégumentaire abondant (poils écail- 
les...) (Entomobryomorphes autres qu'Isotomidae, Symphypléones), 
certaines petites formes de ce deuxième ensemble pouvant bien 
entendu revenir secondairement à un mode de vie interstitiel (Onco- 
poduridae, Neelidae euédaphiques par exemple). 

Le microcosme que constitue le rocher du Lac du Merle traduit 
bien ces deux aspects évolutifs. 

En conclusion, en ne considérant les espèces que sous l’angle de 
leur aptitude à coexister et les relevés sous celui de leur possibilité 
d’héberger des peuplements identiques, on peut dire que deux influ- 
ences principales peuvent rendre compte des affinités : 


— un facteur lié à la nature du milieu (sol, mousse) intégrant 
bien entendu les diverses caractéristiques de ces biotopes. 

— un facteur d'une importance légèrement moindre, en relation 
avec la disponibilité hydrique du milieu, et dont la teneur en eau 
des mousses au moment du prélèvement peut donner une image 
instantanée, reflétant leur histoire hygro-thermique au cours du 
cycle saisonnier. 

Cet exemple de micro-répartition montre le danger qu'il y a en 
écologie édaphique à considérer comme homogènes des stations 
déterminées la plupart du temps sur des critères uniquement phy- 
tosociologiques ou pédologiques, mais groupant en fait des biotopes 
très divers, danger d'autant plus grand qu'on tirera des résultats 
obtenus des généralisations abusives. 

Les micro-arthropodes du sol exigent des techniques d’analyse 
à leur échelle. Quant aux essais de synthèse, ils nous semblent 
pour le moment bien prématurés. 
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